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SUMMARY 


A simplified analysis is made of the shielding mechanism which 
reduces the stagnation-point heat transfer when ablation takes place at 
the surface. The ablation process considered is that in which the mate- 
rial changes directly from the s olid to the gaseous state . The automatic 
~fih-tfi1dTng rn^e hari-T mi Is dismissed and the significant thermal properties 
of a good ablation material are given. The results of the analysis are 
given in terms of dimensionless parameters. 


INTRODUCTION 


It is well-known that the rate of heat transfer to a body is reduced 
when the nose or leading edges are blunt. Such a shape causes a thick 
boundary layer so that a shield of air over the nose or leading edge 
reduces heat transfer to the body. This reduction In heat transfer is 
not always sufficient to achieve desirable body temperatures, however, 
and further shielding Is then required. 

This shielding Is accomplished by providing a material which in some 
way absorbs heat; for example, a copper nose of sufficient thickness 
absorbs enough heat to shield the interior of the body. Alternatively, 
the material may be supplied a6 a gas which is injected through a porous 
wall into the laminar boundary layer near the stagnation point. (See 
ref. 1.) Another method is to allow ablation to take place. This method 
Is not just a negative approach in which heat is expended on the sublima- 
tion of unimportant material; it is also an extension of the blunt-nose 
idea of using the thickened boundary layer to shield the nose. The con- 
vection of heat by the foreign gas Introduced can make the main contribu- 
tion to the shielding at sufficiently high stagnation temperatures. 

Previous theoretical studies of stagnation-point ablation have been 
concerned mainly with melting surfaces in which a layer of liquid flows 
over the surface (refs. 2 and 3)j the main conclusion was that the con- 
vective shielding by the liquid layer Is limited by the temperature at 
which the liquid boils. When vaporization occurs, however, the con- 
vective shielding in the gas boundary layer can increase considerably. 
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Recently, a comprehensive review and study of the convective heat transfer 
with mass addition has been made (ref. 4) which includes the effects of 
dissociation, recombination, and combustion in the boundary layer. A 
reference list on the subject may be found in that report. 

The purpose of the present report is to show that the mechanism of 
shielding by vaporization may be explained quantitatively and qualitatively 
by a simple analysis in which the more important physical aspects of the 
problem are correctly represented. The results of this analysis indicate 
the significant material properties that are required to make a good abla- 
tion shield. 


SYMBOIS 


x coordinate along wall 

y coordinate normal to wall 

z transformed y- coordinate 

Z arbitrary value of z outside boundary layer 

u component of velocity in x-direction 

v component of velocity in y-direction 

U free-stream velocity in x-direction 

C constant in velocity distribution 

T temperature 

T effective mean temperature (eq. (40)) 

W mass concentration of foreign gas 

W effective mean concentration (eq. ( 38 )) 

p density of mixture 

H coefficient of viscosity 

k thermal conductivity 

D-) p coefficient of binary diffusivity 
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L 

S 

°b 

°P 

%u 

R w 

N Sc 

Np r 

q. 

6fji 

% 

0 

K 

Heff 


rate of mass loss per unit area of wall 
latent heat of sublimation 


enthalpy ratio. 


°p,2(^e “ Tv) 

Cp,2(Te " T w) + L + %(T W - T5) 


specific heat of solid 

specific heat at constant pressure 

effective mean specific heat (eq. (^3)) 

Nusselt number 

Reynolds number 

Schmidt number 

Prandtl number 

heat -transfer rate per unit area 
velocity boundary-layer thickness 
thermal boundary-layer thickness 
concentration boundary-layer thickness 
thermal-layer thickness within the nose 
constant defined by equation (39) 

> effective heat capacity of ablation material 


Subscripts : 

e external flow, near stagnation point 

v wall 

b body 


0 


no injection 
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1 foreign gas 

2 air 


SHIELDING MECHANISM 


The configuration considered is that shown in figure 1. The coor- 
dinate system is fixed. in the surface at the stagnation point; thus, the 
interior of the body moves in the y-direction toward the surface with 
velocity v^, the ablation rate. Within the solid in the steady state 

there is a balance of conductive and convective heat transfer. 

At the surface ablation takes place and gases are produced; these 
gases are transported away from the wall and are convected in the laminar 
boundary layer near the stagnation point. This convection of gaseous 
material provides additional shielding and reduces the heat transfer to 
the surface. The rate of mass transfer at the surface is itself deter- 
mined by the rate of heat transfer to the surface; thus this process is 
an automatic shielding mechanism. Additional heat transfer to the sur- 
face causes additional mass transfer with an attendant increase in 
shielding. 

In the steady state the boundary layer will consist of a velocity, 
temperature, and concentration boundary layer, each generally of different 
constant thickness. It is assumed that the component of velocity parallel 
to the wall u is linear in x and that the normal component of veloc- 
ity v, the temperature T, the concentration of foreign gas W, and the . 
properties of the mixture are functions only of y, the distance normal to 
the wall. (These assumptions result directly from the similarity nature 
of the flow and the absence of thermal and concentration gradients along 
the wall. ) 

For a given external flow and stagnation temperature the gas given 
off in the resulting ablation diffuses through the boundary layer and is 
convected with the air, as a mixture, under the action of the pressure 
gradient imposed by the external flow and the shear stress due to the 
presence of the wall. The rate of mass loss and the concentration of 
foreign gas at the wall are uniquely determined by the external-flow 
conditions and the properties of the solid material. 

The overall solution of the problem is obtained by solving separately 
the heat-transfer problem in the solid and the heat and mass-transfer prob- 
lem in the boundary layer; these solutions are then matched by the boundary 
conditions at the ablation surface. 
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AMETSIS 


Heat Transfer in the Solid 

The transfer of heat in the solid moving with constant velocity v-^ 
toward the surface is governed hy the following equation: 




( 1 ) 


Diffusion of heat Convection of heat 

to interior toward surface 


It is assumed that temperature gradients parallel to the surface are 
negligible and that c^ and k^ are constants. Equation (l) holds when 

the thickness of the solid is large compared with the thickness 6 of 
the thermal layer (within the solid and near the surface ) in which the 
temperature changes rapidly from its value at the wall T v (the abla- 
tion temperature) to its value at the far interior 1^. 

The solution of equation (l) which satisfies the conditions at 
y = 0, T = T w (constant) and as y -» -co, T = T^ (constant) Is 


. . k b 

T = Th + (T v - Tb)e 


( 2 ) 


A characteristic thickness 9 is defined by 


e = 



?w - Tfa 



(3) 


thus, the heat content of the body per unit surface area is 0^0(1^ - Tb). 
Thus, in order to confine the high temperatures to a 15a in layer near the 
surface, the material should have low conductivity and high specific 

heat . Equation ( 3 ) also shows that the thermal layer decreases in 
thickness as the rate of mass loss increases. 
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From equation (2) the rate of transfer of heat away from the surface 
at y = 0 is 


k b ~j = °b( T w - *b)m W 

v y /w- 


where w“ denotes that the surface is approached through negative values 
of y. Equation (4) is a statement that the transfer rate at y = 0 is 
that required to raise the mass m of specific heat c^ through the 

temperature difference T w - Tb* The rate of mass loss m is unknown 
at this stage of the analysis and must be determined by consideration of 
the boundary-layer flow over the surface of the wall. 


Wall Conditions for Ablation 
The heat -transfer condition at the wall is 



Im + - Tb)m (5) 


Heat transfer Latent 

to wall heat 


Body 

heat 


The last term of equation (3) is obtained from equation (4) and is the 
heat lost when the temperature of the mass m is raised to tempera- 
ture T w before ablation. 

The mass-transfer condition is 


■n dW 

^ + 

Diffusion of gas 
away from wall 


(3 w v w^ 

Convection of gas 
away from wall 


m (6) 


Gas introduced 
at wall 


where the subscript w refers to conditions that exist as the wall is 
approached through positive values of y. 
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A third condition is that which states that there is no net transfer of 
air into the wall: 




(l - W w )p v v v 


(7) 


Diffusion of air 
toward wall 


Convection of air 
away from wall 


that is, there is a balance of diffusion of air toward the wall because 
of the concentration gradient and convection of air, of concentration 

( dW \ 
p Tr P-| p — I from 

^/w 

equations (6) and ( 7 ) gives Pv v w = nay be taken to be a defini- 

tion of v v . The conditions denoted by equations ( 5 ), (6), and ( 7 ) are 
required for the solution of the boundary- layer equations. 


Boundary-Layer Integral Equations 

The method used in this section to determine the effect on the heat 
transfer of convective shielding in the boundary layer is that formulated 
in reference 1. The principal equations are repeated herein (with some 
modification) for the sake of completeness. 

The boundary-layer equations for two-dimensional flow may be written 
as follows: 


-p v V(Z) + m 



( 8 ) 


Mass flow into 
boundary layer 
from external, 
stream 


Mass flow 
at wall 


Convection of 
mass along wall 


m 



Wgdz 


Gas introduced 


Gas convected 


( 9 ) 
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/ \ z 

-Cp >2 (T e - T„)cvV(Z) - L = C, c p (T - T„)g dz (10) 


Heat flow into 
boundary layer 
from external 
stream 


Heat transfer 
to vail 


Convection of beat in 
boundary layer along vail 


In equations (8), ( 9 ), and (10) C is the constant vhich appears 
in the external-velocity distribution 


U = Cx 


v(z) = — 

Pw 


- r 7 ± 

J 0 p v 


dy 


and Z is some arbitrary value of z outside the boundary layer. 

The quantity V(Z) may be eliminated from equation (10) by use of 
the relation 


C P = °p,l W + " W) 


and equations (8) and (9)j the result is written 


(ID 


CfV f* Pp(T e - T)^ dz = ( k + °P*l (T e " (12) 

Alternatively, by use of the boundary condition of equation (5)> 
equation (12) has the form 

Cp v jT Cp(T e - T)g dz = Jcp,i(Te - T w ) + L + (13) 



MCA TN 4-392 


9 


The equations for transfer of heat and mass near the three-dimensional 
stagnation point also reduce to equations (12) and (13) if x is measured 


from the stagnation point and 


z 



ay 


Heat-Transfer Relations for no Ablation 

The rate of mass loss experienced by the body is determined in part 
by the aerodynamic heat transfer which would prevail in the absence of 
ablation; this rate is denoted by q^. 

The heat-tr an sfer relations here are found in terms of the Nusselt 
number and Reynolds n umb er which have the usual definitions near the 
stagnation point: 


%u,w 



Rv « c £2 X 2 


In addition, the Prandtl number and Schmidt number are defined as 


^Pr,w 


iV 

^w 


C P,2 


K Sc,w 


Pw°12 


The rate of transfer of heat qQ to the wall, when there is no 
ablation, is then given by 


<30 = 





(1M 
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The dimensionless 
erence 5 and used 


quantity 5^ 

R v 1/2 
in reference 1. 


used herein is that suggested in ref- 
Por the axi symmetric flow. 



0.765 


' Pe^e 


k 0.4 

N- 


Pr,w 


0.4 


and for the two-dimensional flow, 



0 . 570 : 


' P e P-e 


v0.4 


^Pr,w ( 


0.4 


(15) 


(16) 


In order to obtain an approximate expression for the thicknessl &u,0 of 

the velocity boundary layer, equation (12) with m = 0 is written in 
the form 




k r 


1 f Z - I u dz . X /M l/2 1 /% u.A 

^u,0 J 0 T e ” ^w u V^PwCy %>r, wL I 


(17) 


and the linear profiles 


u _ z 
U &u,0 




> 

•*~e - T _ 1 z 
T e " ^w ^,0 


(18) 


are used in order to evaluate the integral in equation (17). (The 
quantity 8r ; o I s the thickness of the thermal boundary layer.) 

•4he quantity &u,0 is more correctly related to twice the velocity- 
layer thickness in the axisymmetric case. 
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Equation (17) then has the approximate form (ref. l): 



i /u, N 1 / 2 n Afo \ 

1 / Pw \ 1 /■"Ku,v Y 

0 yPv 0 / N Pr , w\R v l/ 2 y 0 


(19) 


Equations (15) and. (l6) show that the effect of Prandtl number on the 
heat-transfer parameter is given by 


%u,w\ 

<r* • 0 


^%u,v 


(1) 




Npr, 


0.4 


v 


'0 


( 20 ) 


where the superscript (l) denotes the value when Np r ^ w 


= 1. 


When Npj. w = 

3 6m 

similar; thus, — =. 

&u 


1, the thermal and viscous diffusive effects are 
= 1. From equation ( 19 ), 



1/2 . 6 ^ U) 


\Rw 1/2 / 0 


( 21 ) 


is assumed to he independent 


( PyC \^-/2 

The dimensionless thickness 6,, n 

’ V J 

of the Prandtl number; thus, when equations (20) and (2l) are inserted 
into equation ( 19 ) , the following result is obtained: 



= 


-0.5 


In the following section it is assumed also that the relation 


( 22 ) 


°u 


is valid even when ablation takes place. 
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Boundary -Layer Flow With Ablation 

Consider now the effect, on the boundary layer, of ablation at the 
wall. The boundary layer becomes thicker and the profiles may change 
(because of the increase in mass flow) with a resulting reduction of 
gradients in the bo undar y layer. In particular, the heat transfer to 
the wall is reduced since the heat convected parallel to the wall in 
the boundary layer is increased. 

The increase in boundary-layer thickness is found by considering 
the mass flow in the boundary layer in the direction parallel to the 
wall: 


CfVj 

0 u 

= 

r^v.,o u 

°°»Jo S 12 

+ 

m 

Mass 

with 

flow 

ablation 


Mass flow with 
no ablation 


Mass 

ablated 


which gives 




(23) 


when the linear profiles — = — - — 

U &u,0 

ablation, respectively, are used. 


and — = — for no ablation and 
U -B* 


In order to keep the ideas simple, consider first the case in which 
the gas given off has properties identical with those of air. (The more 
general case will then be given as an extension of this method. ) 


The heat-transfer equation (eq. (12)) may be written 


1 : 

r Z - 1 « to - 1 I 

w 'i 

1/2 

I 

1 %u,w m 

&u J 

0 ^e “ ^w U §u' 


KPt,v Ew 1/2 


by use of equation (l4). Equations (l8), (23), and (21) are used to 
simplify this equation with the following result: 
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_1_ . /_1_ .Li % -0.6^ 

HPt,v Ew 1/2 \Xpr,v^p£] o \ 3 Pr ’ V 


m 




(24) 


-I Niff- -rj. 

where is given "by equation (15) or (l6). Equation (24) 

N Pr,v Rv 1 / 2 

rewritten in terms of the heat-transfer rate q^ by using equations (l4) 
and ( 15 ) and the relation 


q. 


c p ,2(T e - T w )( lVJv C) 1 / 2 


1 -%UjW 

N Pr,w Rw l/ 2 


is 


qo = 


°t>(^w “ + L + Cp^2(^e “ ®w)^l “ ^ ^Rr,w 


m (25) 


Body 

heat 


Latent 

heat 


Boundary-layer 

shielding 


The quantity 


1 - 0 . 6 1 / , 

1 - •=• J(T e - T v ) is the effective, temperature 


rise in the boundary layer of the gas given off in ablation. 

When qQ is expressed in terms of dimensionless quantities by use 
of equations (l4) and (20 ) , equation ( 25 ) becomes 


c p,2( T e - T-w) 


Npr, 


0.6 


m 


w 




1/2 


Cp^2(^e ~ Tw) + L + c^/Ty- - Tfc,) 


^%u,w 




( 1 ) 

r 0 


( 26 ) 
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Equation (26) may be used to determine m . ■ when the enthalpy 


ratio 


c p,2(l'e “ T w ) 

Cp^C^e ~ *w> + L + o^(T w - T-jj) 


is known. The result is 




\Nu,w 


^ l/2 /o 


, T -0 . 6 _ 
Npr,w- S 

T 1 -NT “ 0 . 6 _ 

1 _ _ Np r ^ w s 


As a check on the validity of the several simplifying assumptions 
that have been made, a comparison is made with exact solutions for air- 
to-air injection which were compared with the present method in ref- 
erence 1 for axi symmetric flow. When ablation occurs, the amount of 
gas introduced into the boundary layer is no longer a free parameter, 
as in the case of injection, but is determined by the rate of heat 
transfer according to the following equation: 


1 %u,w _ L + °b( T w ~ T b) m 

Npr,w R w l/ 2 Cp, 2 (Te - T v ) 


( 28 ) 


which is the dimensionless form of equation (!?)• 


Thus, the dimensionless rate of mass loss 


by the intersection of the straight line 


is determined 


(Pw+V^) 


%u,w _ 1 - S m 

„ 1/2 Pr ' v S , .1/2 

(fW c ) 
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(which is eq. (28) by definition of S) with the straight line given by 
equation (24) . This intersection has been determined for various values 
of S and the results are shown in figure 2. 

The agreement between the results of the present analysis and the 
exact air-to-air injection results of reference 6 is extremely good 
except for values of the parameter S very near unity. The latter 
situation occurs when the latent heat L and the body heat %(T W - Tq) 

are small. Such a situation is undesirable; thus this range of S is 
unimportant from a practical viewpoint. The foregoing comparison iB 
considered to be sufficient justification to extend the method to take 
account of the additional effects which occur when the gas has properties 
different from those of air. 


The mechanism by which air diffuses through the boundary layer in 
a binary mixture is exactly analogous to that by which heat diffuses in 
the absence of viscous dissipation when Cp^q. = Cp^g* The analogy is 

seen when the gas -transfer equation ( 9 ) is written 




m 



(29) 


and compared with the heat-transfer equation when Cp^ q = Cp^ 2 



(50) 


obtained by the use of equation (13) and the definition of S. 

The boundary condition at the wall which is obtained from equation ( 7 ) 




1 - w, 


w 


W, 


m 


w 


w 


(31) 
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is compared with 


P 


k d foe ~ T 
pe p ,2 <3y\T e - 


v 



( 32 ) 


which is obtained from equation ( 5 ) and the definition of S. The 
analogous quantities are given in the following table: 


Binary diffusion 

. Thermal diffusion 

W 

c p,2(Te “ T) 

Ob(Tv - Tb) + L + Cp^2(Te - T w ) 

d 12 

k 

PCp,2 

Hsc " p£ 2 

Npr = ^ cp,2 

W w 

c p,2(Te - T v ) _ g 

Oja^Tv ” Tb) + L + Cp^2(Te “ T w ) 


Because of this analogy the following relations may be written: 
From equation (22), 


6W 

^u 



( 33 ) 


where is the thickness of the concentration boundary layer, and 

from equation ( 2 6 ) 
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N Sc,w 


0.6 


m 


W v = 


(Pw^C ) 1 / 2 


Nlfu,w 


,(D 




1/2 


+ 1 

3 


m 


(Pw^ c ) 


1/2 


(34) 


where 


%u,w 
p V 2 



is given by equation ( 15 ) or (l6). 


Equation (34) gives the concentration of foreign gas at the wall, 
which is required for the determination of the boundary-layer shielding 
(since the specific heat of the mixture depends on W). The heat-transfer 
equation (12) is rewritten by using equations (5) and (ll) as 



and the linear profile 


¥ 

Ww 



(36) 


is assumed. 

Substitution of equations (l8) and (36) into equation (35) gives 
an equation analogous to equation ( 25 ) but includes the effect of the 
difference in specific heats and a binary diffusion (in the quantity 
N Sc,w) • This result is written in the form (see ref. 1 for a more com- 
plete derivation) : 
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«o 


Qb^ " ®b) + 


+ 


Body- 

heat 


latent 

heat 


|°p,l^ + {^- " J H Pr,w j( T e “ T w) 


Boundary-layer 

shielding 


(57) 


where W is an effective concentration of foreign gas in the boundary 
layer and is defined by 


W = 


1 " *%c,v 


0.6 


1 " 5 


- 0.6 


( 58 ) 


where 


K 




JL T e ~ J 

W w T e - T w 



( 39 ) 


Again the effective temperature 
layer is ^ “ 5 N Pr ,w"°' ^( T e “ 
convected gas is 


rise of the convected gas in the boundary 
T w ) and the effective temperature of the 


T = T w + ^1 - | Npr"°’ 6 j(Te - T w ) (^0) 

which increases with Np^y.. (See fig. 3.) 

When is expressed in terms of the enthalpy difference across 

the boundary layer by use of equations (l^-) and (15) or (l6 ) , the fol- 
lowing expressions are obtained for the rate of mass loss: 

For axlsynmetric flow. 
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°^ 65 fe) N Pr ) v'° -6 (lW c ) l/ %,2( T e - %,) 


°b(Tw - Tb) + 1 + Cp(T - GV) 


and for two-dimensional flow, 


0.570 


PeM“' w -0.6/ n .l/2 

q. jx ] N Pr/W (pv^V-O °p,2( T e 


Pw^v , 


°b(lw - H?b) + L + Cp(T - T w ) 


where 


C P = C PA W + °p,2(! " w ) 


is the effective specific heat. 

Equations (40) and (ll) show that 


°* 765 fe) (PwPwC) l/2 WPr,w"°- 6 


°p L 1 w - 0.6 

STil " 5 Pr ’ w 


for large values of T e - T v . The effective total heat capacity of the 
shielding material is given by equation (37) as 


H eff - = %(T V - Tb) + L + 3p(T - T w ) 


which increases with Te “ T w . 
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DISCUSSION 


It is seen from equations (4l) and (40) that the rate of mass loss 
is reduced by choosing a material such that the enthalpy ratio 

OfaC^W - Tb) + L + Cp(T e - Tw)(l - | Npr,v’°* 6 ) 


Cp,2(^e " ^v)Npr^v 


- 0.6 


is as large as possible; this condition requires that the following 
quantities be large: 

(a) The surface ablation temperature T w 

(b) The specific heat of the solid 

(c) The latent heat L 

(d) The effective specific heat Cp 

(e) The Prandtl number Np^^. 

In addition, as shown in equation (3) , the thermal conductivity k^ of 

the solid material should be small in order to confine the high surface 
temperature to a thin thermal layer near the surface. 

The effective specific heat Cp is large when > °p,2 811,1 w 

is large. Since W increases as %c,w decreases and is virtually 
independent of Np^^. for 0.6 < Npr,w < 1*6 (fig. 4), "Nsc,w should 
be as small as possible. When Cp,i < Cp,2> however, W should be small 
and hence Ns C;W should be large in order, to reduce the rate of mass 
loss. These diffusive effects on Cp are explained physically in ref- 
erence 1 and may be summarized briefly as follows: When °p,l > °p,2 

and Nsc,w is small, the foreign gas diffuses quickly through the boundary 
layer before being converted; when Cp ; i < °p,2 811(1 % c ,w ls lar se^ 

the gas remains near the wall but displaces the air away from the wall. 

In both cases the fluid of greater specific heat is convected in that 
part of the flow at the highest temperature and with the highest velocity, 
and thus produces the most convective shielding. (See fig. 5«) 
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, C P,2 


In order to bring out the more important shielding parameters 

and without reference to any particular material, the Prandtl 


number and Schmidt number are given the value of 1; whereas in a more 
general calculation they would be determined by the concentration W 
and the properties of the material in the gaseous state. (See ref. 1.) 
A further simplification is brought about by the approximation 


Iwv) C( W C)1/2 “ (»e^C ) l/2 


in equations (4l) and (42). Figure 6 shows the rates of mass loss as a 


function of 


c p ,2 


and S for axisynmetric and two-dimensional flow; 


these rates were determined by using these assumptions. 


CONCLUDING REMARKS 


The cooling of a blunt body by ablation appears to be a promising 
method since the shielding mechanism is such that gas shield increases 
as the heat transfer to the body increases. Furthermore, the effective 
heat capacity of material increases even though the rate of mass loss 
becomes almost constant as the heating rate increases. The most 
important thermal properties for a good ablation material are: low 

thermal conductivity, high specific heat of solid, and high specific 
heat in the gaseous 6tate. The simple quantitative results of the 
analysis bear out the physical ideas of the ablation mechanism and 
agree well with the results of an exact analysis for an air-to-air 
injection. 


Langley Aeronautical Laboratory, 

National Advisory Conmittee for Aeronautics, 
Langley Field, Va., July 3> 1958* 
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Figure 2.- Rate of ablation when foreign gas is same as air. 3?wo- 
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Figure 3»- Effect of Npj. v on T 
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Figure 5-- Schematic diagram of diffusive effects. 
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(b) Two -dimensional flow; Np*. v = 
Figure 6.- Concluded. 
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